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Thermal right-handed neutrino self-energy 
in the non-relativist ic regime* 

M. Laine 

ITP, AEC, University of Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland 

Recently the issue of radiative corrections to leptogencsis has been raised. Considering 
the "strong washout" regime, in which OPE-tcchniqucs permit to streamline the setup, we 
report the thermal self-energy matrix of heavy right-handed neutrinos at NLO (resummed 
2-loop level) in Standard Model couplings. The renormalized expression describes flavour 



O j transitions and "inclusive" decays of chemically decoupled right-handed neutrinos. Although 



CP-violation is not addressed, the result may find use in existing leptogenesis frameworks. 



1. Introduction 



Q\ , Leptogenesis is currently among the most popular scenarios for explaining the observed cos- 



mological baryon asymmetry. Surprisingly, although the basic mechanism is fairly simple [1] 



OO; 

it appears difficult to develop a fully consistent theoretical description of the physics involved. 
The reason is that many different subtle topics, such as CP-violation, baryon number viola- 
| tion, deviations from thermal equilibrium, as well as resummations necessary for a systematic 

treatment of relativistic thermal field theory, need all to be put consistently together. 

As far as CP-violation is concerned, it could be both of "direct" and "indirect" type, 
^ ' known as "vertex" and "wave function" corrections, respectively [2, 3, 4]. The indirect 

CP-violation originating from flavour oscillations is non-trivial even at low temperatures 
where essentially the vacuum formalism can be used [5, 6]. Moreover, collective thermal 
phenomena and higher-order corrections could be important [7, 8, 9]. At present many 
approaches remain phenomenological, but efforts towards a more systematic treatment are 
under way [10]-[14].t Ultimately the goal should be to present theoretically consistent results 
in terms of Standard Model couplings, in the sense that have recently been obtained for 
CP-conserving rates both in the "ultrarelativistic" (m to p ^ M <C irT) [18, 19] and "non- 
relativistic" (m top < ttT <C M) [20, 21] regimes. Here we concentrate on the latter regime, 
called the "strong washout" case in that memory about initial conditions has been lost. 



* Presented at the International Workshop Frontiers in Perturbative Quantum Field Theory, 10-12 Septem- 
ber 2012, Bielefeld, Germany. 

^The situation is analogous to that for CP-violation in phase transition-based baryogenesis [15, 16, 17]. 
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In efforts towards systematic leptogenesis, the right-handed neutrino self-energy plays an 
important role [22, 23, 24]. Examples of recent discussions can be found in sec. 3 of ref. [13] 
and in sec. 4.1 of ref. [14], in both of which the self-energy was handled at leading order in 
Standard Model couplings. In the temperature regime of interest, right-handed neutrinos are 
out of equilibrium, but the Standard Model particles are in equilibrium, at least as far as 
CP-conserving reactions are concerned. Therefore the self-energy of the right-handed neu- 
trinos, which reflects the dynamics of the other particles, can be computed with established 
techniques of thermal field theory, and there is no reason to restrict to leading order. 

In this note a result for the thermal right-handed neutrino self-energy in the non-relativistic 
regime is presented at NLO (partly even NNLO) in Standard Model couplings. Compared 
with our earlier work [21], where the production rate of right-handed neutrinos was computed, 
the real part of the self-energy is added here, and the full Lorentz and flavour structures 
are included. Technically, we work at order 0{h^ u h u ) in neutrino Yukawa couplings, thus 
addressing CP-conserving processes, whereas CP-violation originates at the order 0(hlh u ) 2 . 
As an outlook, the same observable could in principle also be computed in the "relativistic" 
(irT ~ M) and ultrarelativistic (irT 3> M) regimes; its physics is related to "CP even 
damping and time evolution", as discussed e.g. in sees. 5.2-3 of ref. [25]. 



2. Self-energy matrix 

We start by computing the Euclidean correlator (the Lagrangian and other conventions are 
summarized in appendix A) 

V B (K) = J P dr Jj K - x ^aJ){X) (£aJ)(Q)) T , (2.1) 

where X = (r, x); K = (Av^k), where k n are fermionic Matsubara frequencies; t is a lepton 
doublet; and (j> = io2<t>* is a Higgs doublet. The corresponding retarded correlator is obtained 
through the analytic continuation k n — > —i[k° + i0 + ]. 

Treating 7 5 in Naive Dimensional Regularization and employing Feynman gauge (other 
regularizations and gauges were discussed in ref. [21]), we obtain the contributions (dashed 
line denotes Higgs, solid lepton, double top, wiggly VK ± ,Z°) 
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p 



P 2 (K-P) 2 ' 



(2.2) 




r PQ Q 2 P 4 (K-p) 2 ' ( } 



2 



2 , o„2 ^ D ^ ia L ($T - f )a F 



- (Sib + 3 #2b) 2" ^ Q2p4(K - P) 2 ' '' 2 ' 5 ' ) 
<V> = (9 lB +3g 2B )-^ piQ2{p _ Q)2{K _ p)2 , (2.6) 
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Here D = 4 — 2e; ^ Q are bosonic sum-integrals; and is a fermionic one. 

The sum-integrals in eqs. (2.2)— (2.8) are of a "tensor" type [26]; the numerator transforms 
non-trivially in 0(4) rotations. It is known that in an ultrarelativistic plasma, fermion 
self-energy does not respect Lorentz invariance; the part multiplying 7 differs from that 
multiplying the spatial -y k [27]. It turns out, however, that in the non-relativistic regime, 
i.e. (vrT) 2 <C K 2 , in which the problem can be discussed with OPE language [28], the first 
thermal corrections do respect Lorentz invariance. Perhaps the simplest way to understand 
this is that thermal corrections of 0(ttT) 2 are represented by the condensate {4>^<j)) T [21], and 
there is no room for tensor structures in this condensate. 

In practice, the sum-integrals are of two types. Those proportional to the external four- 
momentum $T are scalars; results can be found in appendix C of ref. [21]. Those proportional 
to ]/> , Q or Jji are tensors; Matsubara sums need to be carried out separately for the temporal 
and spatial components. In the OPE regime this can be done with the techniques introduced 
in refs. [29, 30]. In each case, only structures proportional to $T survive at 0(7rT) 2 . 

Concretely, multiplying Tie by the renormalization factor related to the neutrino Yukawa 
couplings, Z v = 1 + ^y2 e [\ht\ 2 N c — § (<7 2 + 3g 2 )] + 0(g 4 ) where g denotes generic Standard 
Model couplings and corrections involving neutrino Yukawa couplings were omitted, we obtain 

Z„E E {K) = a L i^ a B I -^(- + ln^ + 2 
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Here Z m = 1 + u^yri [6A + \h t \ 2 N C — | (gf + Sgf)] + C(5' 4 ) is the renormalization factor related 
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to the Higgs mass parameter (denoted below by to 2 ,) and, for 7rT>m t0 p [21], 

rp2 
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+ 0(g 3 ), 
(2.10) 

(2-11) 



Here m D1 , m D2 are the U Y (1) and SU L (2) Debye masses, and chemical potentials have been 
set to zero. All thermal corrections lie in Z m (<f>' 4>) T . If we count (nT) 2 ~ K 2 then eq. (2.9) 
is "NNLO" because, after analytic continuation, its real part contains terms of O(g ), 0(g 1 ), 
0(g 2 ) (cf. eqs. (2.10), (2.13)), and its imaginary part terms of O(g ), 0(g 2 ), 0(g 3 ) (cf. 
eqs. (2.10), (2.14)). However, in the non-relativistic regime it is more natural to count 
(irT) 2 ~ g 2 K 2 , and then the imaginary part is complete only up to NLO. 

For the next steps, a formalism needs to be chosen by which to represent right-handed 
neutrinos; various possibilities are summarized in appendix A. If the neutrinos are represented 
as chiral Dirac fermions then, after the analytic continuation k n — > —i[k° + i0 + ], eq. (2.9) 
amounts to a "wave function correction" in their retarded self-energy: 

a L ^a R ^a L ^a R {l + / l t/i iy [0 R (/C 2 ) + isign(fe o )0 I (/C 2 )]} , (2.12) 

where /C = (A;°,k), h v is a 3 x 3 matrix of Yukawa couplings (cf. eq. (A. 5)), and 
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10) 



(2.14) 



(4vr) 2 K 2 

If we rather employ Majorana spinors, then eq. (A. 7) implies that the flavour structure can 
be "reflected" to the left-handed components: 

Yf. {l + [^t^ + ^f a L ][</. R (/C 2 )+zsign(A ; o )0 I (/C 2 )]} . (2.15) 
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Given that h\,h v is Hermitean, we denote (h\,h u ) T by (h\,h v )* in the following. If a time- 
ordered correlator (X T ) is considered rather than a retarded one (Er), then the corresponding 
self-energy reads 



where ra F denotes the Fermi distribution, and Im refers to a cut across the k -axis. In the 
non-relativistic regime, \k°\ 3> ttT, we can simplify 1 — 2n F (/c°) = sign(/c°), and therefore in 
the time-ordered case it is the combination <p n (K?) + ^(/C 2 ) that appears. 

3. On-shell renormalization 

As suggested by eqs. (2.12), (2.15), it is conventional to view the divergences appearing in R , 
cq. (2.13), as being cancelled by wave function renormalization. Wave function normalization 
being somewhat unphysical, it may be more elegant to discuss only mass renormalization 
explicitly. This can be achieved by considering renormalization at an on-shell point. 

In the literature, various formalisms are being used for addressing right-handed neutrino 
dynamics in the context of leptogenesis. Within the "Kadanoff-Baym" framework, retarded 
and advanced self-energies appear directly, and the pole positions of the corresponding prop- 
agators can be solved for (cf. e.g. sec. 4.1 of ref. [14]). Within the "canonical" approach, a 
density matrix between free one-particle states is constructed and evolved [31, 32, 33]. In 
both cases on-shell particles appear and the considerations below may apply. 

Of course, strictly speaking right-handed neutrinos are not on-shell states, because they 
decay. Nevertheless, the width being parametrically suppressed by 0(hlh u ), the concept 
of an on-shell particle may still be a useful practical notion, like for the top quark. In the 
following, we refer to the pole mass as the real part of the pole position in the A;°-plane. 

Let M be the renormalized mass matrix which we choose to be non-negative and di- 
agonal, and write M B = M + SM . The diagonal elements of M are denoted by mf. 
M = diag(mi, 1712,1713). Making use of Majorana spinors, the effective action for time-ordered 
correlators has the form 



E T (/C) = Re E fl (/C) + i [l - 2n F (fc )] Im £ fl (/C) , 



(2.16) 




+ [4l lA + (/ i tyX]^[^(-9 2 ) + ^(-9 2 )]}%) 1 *=(t,x), (3.1) 



where SM denotes the mass counterterm. Within the 0(hlh u ) part we can make use of 
tree-level equations of motion; denoting 



M = diag(mi, m 2 ,m 3 ) , rhi = mi 



)] 



(3.2) 
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and recalling from eq. (A. 6) that mass matrices can always be symmetrized, we get 
S cS = J ' -N{X){ip -M-5Ma n -5M*a^ 

+1 [ h l K M + M{hlh v )*] a R + \ [(htK)*M + MhtK] a L }N(X) . (3.3) 

In the on-shell scheme 5M is chosen to cancel the vacuum part of the 0{fv v h v ) correction; 
hence, writing <f> R = (f> R Y) + (f> R , where 4> R corresponds to the three first lines of eq. (2.13), 

SM nn = \ [{hlK) ili2 m i2 ^(ml) + (h\h v ) t2ll (m 2 h )] . (3.4) 

The finite remainder, which we parametrize through 

M r . = diag(m ril ,m r)2 ,m ri3 ) , m ri = m, [^(mf) + i&(mf)] , (3.5) 

contains thermal "mass corrections" (from (jy R ) as well as a non-Hermitean part, representing 
vacuum and thermal "decay widths" (from : ). 

To be concrete, consider the canonical formalism and define a free Hamiltonian as 

#S (*) = \ ■ HX) {~il k d k +m)n{X): . (3.6) 
From eq. (3.3), the first-order correction reads 

- |r£(*)= -\ : N(X)(^{Re(hlh v ),M r } + % ^[Tm(hth v ),M r ])N(X) : . 

(3.7) 

Defining a state with a specific flavour (ii), momentum (ki), and spin (n) formally as 

|nk 1 r 1 )^(2vr) 3 / 2 al ikiri |0), (3.8) 

and inserting the interaction picture field operators from eqs. (A. 10), (A. 11), matrix elements 
can be computed. Standard properties of the on-shell spinors (cf. appendix A) lead to 

(iiknl #2(0) \i 2 kr 2 ) = S ili2 S TlT2 /C° , /C° = Jtf + ml , (3.9) 



amounting to a free one-particle energy. For i\ = i 2 , the 0{h\,h u ) correction yields 



,2 



(iiknl W«(0) Kikr 2 ) = -^LR^tfcJ^ <5 rm ^ T) «) , (3-10) 



n 

which may be interpreted as a thermal mass correction: <5/C° ~ Re(/i^/ij / )j 1 j 1 2 m (^^) T //C^ ) 1 . 
The decay width is 

<iikri|f2(0)|iik75) = -p^Re^/iJ^V^iK)' ( 3 - n ) 

agreeing with refs. [20, 21]. Flavour non-diagonal parts are non-zero and can be worked out 
from eqs. (A. 16), (A. 17); in the non-relativistic regime, k <C rrii, the spin-part coming from 
7 5 is seen to be suppressed, and the flavour structure is given by the first term of eq. (3.7). 
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4. Summary 

Addressing the full problem of leptogenesis is demanding: although non-perturbative equa- 
tions have been written down, no general systematic solution is known even to leading order 
in Standard Model couplings. However, simpler subproblems can be studied in a controlled 
fashion, and hopefully a similar progress ultimately permeates the whole topic. 

In this note, one simple subproblem has been considered, that of 0{h\,h v ) flavour tran- 
sitions and inclusive decays of heavy right-handed neutrinos in the non-relativistic regime 
( m top^ 7T T <C M). In particular, the NLO corrections to the right-handed neutrino self- 
energy matrix have been determined, cf. eqs. (2.13)— (2.15). Technically, this requires consid- 
ering "tensor-type" sum- integrals [26]. Conceivably, generalizing techniques from ref. [34], a 
similar computation can in the end also be carried out in the relativistic regime. 

The main physics conclusion of this study is that in the non-relativistic regime thermal 
effects are (only) power-suppressed, but nevertheless infrared safe up to the order studied, 
both as far as mass corrections and decay widths are concerned. Whether this statement 
also applies to the actual lepton asymmetry generation is, however, unclear at present. A 
potential obstacle is that OPE techniques [28] , which help to understand the structure of the 
current results as well as the nature of higher-order corrections, are not trivially applicable. 
The reason is that, particularly in the resonant case, CP-violating observables contain the 
heavy mass scale M also in internal propagators, not only in external states. 

Acknowledgements 

I am grateful to Y. Schroder for collaboration on ref. [21] and for helpful discussions. 
Appendix A. Conventions for Majorana spinors 

For completeness, we specify in this appendix our conventions for right-handed neutrinos. As 
usual, they can be represented by two-component Weyl spinors (£), four-component Majorana 
spinors (N), or right-handed chiral projections of Dirac spinors (z/ R = a^y). 
Choosing the Weyl representation for the Dirac matrices, 



where a k are the Pauli matrices, the chiral projectors a LjR = (1 =F 7 5 )/2 become diagonal, 
and the four-component Dirac spinor can be decomposed as 




(A.l) 




(A.2) 
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The charge conjugation matrix can be defined to be 

C ee <7 V= I I ■ (A-3) 





and if we set 

iV - ( -f ) , N={$- e T - 2 ) , (A.4) 

then it is easy to see that N c = CN T = N, i.e. that TV is a Majorana spinor. 

In terms of Dirac spinors, the bare Lagrangian of the right-handed neutrinos reads 

C ee ujj^d^ - (lajj> h UB u R + - u c R M B u n + H.c.) . (A.5) 

Issues related to 7 s were briefly recalled in ref. [21]. Both h VB and M B are matrices in flavour 
space; the Grassmann nature of the fields and the antisymmetry of C imply that the matrix 
M B is symmetric, M B = M B . Making use of Weyl spinors, the Lagrangian can be rewritten 

as 

C = guPd^ - (it 4>h UB £-±e io 2 M B e + H.c.) , (A.6) 

where = (l,crj.), and £ L is to be interpreted as a 2-component spinor. Finally, for the 
Majorana representation, eq. (A.4) implies that 

Ni-fd^N = guPd^, Ni-fd^N = fia^d^*, (A.7) 
NM B a R N = -i T ia 2 M B i, NM*a L N = ^ia 2 M*C- (A.8) 

Omitting a total derivative from the kinetic term (or symmetrizing that in eq. (A.6)), 

£ = I N{i^d^ - M B a R - M B *a L )iV - (la J K B N + H.c.) . (A.9) 

A pleasant feature of the Majorana formulation is that the free on-shell condition has the 
form of the usual Dirac equation, {i^d^ — M B a R — M B a h )N = 0. If the mass matrix is 
diagonal, with eigenvalues rrii G R + , then field operators of the interaction picture can be 
expanded as 

Ni{X) = [ d ' k ^(w^+ BikT 4e^), (A.10) 



where v ikr = CuJ kr and v ikr = uJ kT C; r enumerates the spin states; and K,® = y fc 2 + m 2 . 
The creation and annihilation operators are assumed to satisfy 

{«u kl n4 2 k 2T J = <U<W (3) (ki " k 2) > (A-12) 
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and the on-shell spinors obey — rrii)u ikT = (/£i + mi)v ikT = 0. Writing 

u ikr = f' + m ' Vt , V T = 4f ( ' T M , (A.13) 
^ + m, ^2 V \r) J 

where |r) are eigenstates of n • a with n an arbitrary unit vector, on-shell spinors obey the 
same relations as in the Dirac case: 

~~ VikT2 V ilrn = v-iV.Ti u i\iT2 = ^rrii S TlT2 , u ikri v ikT2 = v ikT1 u ikT2 = , (A. 14) 

"ikn 7°«ikTa = ^kri7°^k^ = 2/C °<W > «ilwi'Ai-l«2 = ^kri7°^-k^ = ■ ( A - 15 ) 

Making use of ^2 T i] T rj T = \ (1 + 7°) 5 the usual completeness relations are readily verified: 
St ^ikr^ikT = + j St ^ikr^ikT = ^i~^- m i- F° r different flavours, the representation 
of eq. (A.13) can be used for showing that 

U UkTl n i 2 kT 2 - T1T2 \ 7= I =— ( > (A.10J 



= H- k NjV^T^--V^T^-^ (A - 17) 
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